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1.  Introduction 


This  technical  memorandum  contains  the  results  of  matched  field  processing  on 
200  Hz  CW  tone  propagation  data  collected  during  the  September  1987  Vertical  Line 
Array  (VLA)  experiment  in  the  North-East  Pacific.  The  Marine  Physical  Laboratory 
(MPL)  digital  array,  long  of  900  m  and  filled  with  120  hydrophones  with  7.5  m  interele- 
ment  spacing,  was  deployed  for  several  weeks  trcm  the  R/P  FLIP,  approximately  400  km 
southwest  of  Monterey,  California.  The  R/P  FLIP  was  in  a  three-point  moor  at  34°58.69' 
N  125°58.43'  W.  Previous  work  [Tranl989]  dealt  with  the  conventional  processing  of  the 
200  Hz  CW  tone  propagation  data  recorded  by  the  VLA,  first,  when  the  source  ship  kept 
a  fixed  stations  89  nm  from  the  R/P  FLIP  on  Julian  271  and,  second,  when  the  ship  towed 
the  source  away  from  the  R/P  FLIP  for  ten  hours  on  Julian  272.  In  both  cases,  the  top  of 
the  VLA  was  at  a  nominal  depth  of  400  m.  The  acoustic  source  was  an  HX90  deployed 
from  the  USNS  NARRAGANSETT.  It  projected  the  200  Hz  CW  tone  with  a  nominal 
source  level  of  184  dB  ref  1  uPa  at  1  m. 

The  results  of  the  conventional  processing  and  of  the  acoustic  modeling  of  the 
propagation  presented  in  [Tran  1989]  demonstrate  the  good  quality  of  the  data  and  the 
ability  to  model  the  propagation  with  the  available  environmental  information.  Based  on 
these  results,  it  is  natural  to  attempt  matched  field  processing  on  the  data  set.  Since  the 
source  level  is  high  in  this  propagation  experiment,  detectability  is  not  an  issue  here.  The 
signal  is  easily  detected  on  individual  phones  using  conventional  techniques.  Therefore 
the  focus  is  on  the  range  and  depth  localization  capability  of  matched  field  processing. 
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2.  Outline  of  the  Matched  Field  Processing 

Conventional  beamforming  can  be  interpreted  as  a  spatial  matched  filtering 
operation.  The  Bartlett  processor,  for  a  particular  frequency  /  is  given  by 

Pb=-J-EhRE  (2.1) 

M 

where  R  is  the  array  covariance  matrix  at  the  frequency  /  based  on  the  data  given  by 

R  =  ^Vx,X,w  (2.2) 

where  X,  is  the  vector  of  Fourier  coefficients  across  the  array  at  the  selected  frequency  / 
corresponding  to  the  time  snapshot  i.  In  Equation  2.1,  E(9)  is  a  normalized  replica  vector 
containing  the  phase  variations  across  the  array,  assuming  a  given  arrival  angle  and  a 
wavefront  shape.  The  replica  vector  E  is  such  that  Ew  E  =  l.  The  beamformer  searches 
peaks  in  angle  to  find  arrivals.  Such  processor  can  be  generalized  to  search  spaces  other 
than  angle.  The  matched  field  processing  techniques  search  in  the  range-depth  space,  that 
is  the  replica  vector  E  is  now  E (r,z)  [Baggeroerl988],  Matched  field  processing  consists 
in  measuring  the  match  or  the  mismatch  of  a  field  probed  by  a  vertical  array  with  a 
predicted  field  for  a  particular  source  range  and  depth.  Matched  field  processing  provides 
detection  as  well  as  localization.  The  main  difficulty  is  to  compute  a  realistic  prediction 
of  the  pressure  field  for  a  source  and  receiver  geometry.  The  replica  vectors  are  com¬ 
puted  by  propagation  models  such  as  normal  mode  models  where  environmental  infor¬ 
mation  is  critical.  The  literature  provides  multiple  examples  of  the  great  sensitivity  of 
matched  field  processing  to  an  imperfect  knowledge  of  the  environment  parameters  such 
as  sound  speed  profile,  geoacoustic  parameters  of  the  bottom,  water  depth  [Porterl987, 
Del.Bazol988,  Feuilladel989,  Tolstoy  1989,  Baxleyl989].  The  acoustic  modeling  per¬ 
formed  in  [Tranl989]  is  used  with  the  ATLAS  normal  mode  model  [Gordonl984]  to 
create  the  replica  vectors  using  a  sound  speed  profile  based  on  the  CTD  and  the  XBT 
data  collected  at  the  R/P  FLIP  during  the  experiment.  The  range-depth  search  space  is 
sampled  every  250  m  in  range  and  every  5  m  in  depth  up  to  a  maximum  depth  of  600  m. 
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Several  processors  can  be  used  in  matched  field  processing  [Fizelll987].  In  addi¬ 
tion  to  the  conventional  Bartlett  processor  given  in  Equation  2.1,  with  E(r,z),  the 
Minimum  Variance  (MV)  processor  [Caponl969],  which  yield  more  resolution  but  is  less 
robust,  is  used  in  the  following.  The  MV  processor  is  given  by 


l  l 


^■SFFi  <2'3) 

where  R  is  defined  like  in  Equation  2.1,  and  E  is  the  replica  vector  corresponding  to  a 
particular  source  range  and  depth.  E(r,z)  is  normalized  to  unit  norm.  In  order  to  be 
inverted,  the  array  covariance  matrix  R  must  be  stabilized,  that  is,  its  main  diagonal 
loaded  by  a  small  amount  of  spatially  white  noise  [Caponl969].  The  matrix  R  is  aug- 

mented  by  y— ^  --I,  where  i  is  the  identity  matrix,  tr  the  trace  operation,  M  the  number  of 


sensors  and  y  the  fraction  of  white  noise,  also  called  stabilization  factor  with  values 
between  10~3  and  KT1.  The  selected  values  for  the  stabilization  in  this  study  are  1(T3,  l(T2 
and  10"1. 


A  preliminary  step  is  to  estimate  the  array  covariance  matrix  at  the  transmission 
frequency,  200  Hz.  The  covariance  matrix  is  computed  by  averaging  126  dyad  products 
X,  Xj1  in  Equation  2.2.  The  array  Fourier  vectors  X,  are  derived  from  two  blocks  of  64  K 
data  points  (1  K  =  1024  point)  using  50  %  overlapped  2048-point  FFTs  (244  mHz 
binwidth),  and  a  Kaiser-Bessel  window  with  a  parameter  of  2.5  (first  side-lobe  level  at 
-57  dB).  Each  covariance  matrix  corresponds  to  4  minutes  and  22  seconds  worth  of  data. 
The  2048-point  FFT  corresponds  to  4  seconds  worth  of  data  and  yields  a  positive  signal- 
to-noise  ratio  on  an  individual  hydrophone. 
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3.  Fixed  Station  Matched  Field  Processing 

3.1.  Description  of  the  Experiment 

On  Julian  Day  270,  the  200  Hz  CW  tone  was  projected  for  various  time  lengths 
and  various  source  depths,  89  nm  (or  165.5  km)  from  the  R/P  FLIP.  The  first  transmis¬ 
sion  with  the  source  at  a  nominal  depth  of  300  m  began  at  20:50  GMT  and  lasted  for  15 
minutes.  The  corresponding  data  were  recorded  on  Tape  915  (which  starts  at  20:48:04 
GMT).  The  source  was  raised  to  a  depth  of  150  m  and  the  transmission  resumed  at  21:30 
GMT  for  10  minuies.  The  data  were  recorded  on  Tape  917  (which  starts  at  21:34:59 
GMT).  The  source  finally  was  raised  to  a  depth  of  20  m  and  projected  for  18  minutes 
from  21:57  GMT.  The  data  were  recorded  on  Tape  918  (which  starts  at  21:58:27  GMT). 

Environmental  data  were  collected  on  the  R/P  FLIP.  Swell  heights  on  Julian 
Day  270  were  visually  measured  between  5  and  10  feet.  A  northern  wind  was  blowing 
with  speeds  between  20  and  25  kts.  XBT  casts  were  collected  at  the  R/P  FLIP  before  and 
after  the  transmissions 

3.2.  Fixed  Station  Modeling 

The  modeling  of  the  propagation  data  was  described  in  [Tran  1989]  for  the  three 
different  depths  at  which  transmissions  were  performed  using  both  the  Generic  Sonar 
Model  (GSM)  [Weinberg  1985]  and  the  ATLAS  normal  mode  model  [Gordonl984], 
Information  relevant  to  the  calculation  of  the  replica  vectors  with  the  ATLAS  model  is 
presented  here,  more  detailed  information  can  be  found  in  [Tran  1989]. 

3.2.1.  Sound  Speed  Profiles 

As  mentioned  earlier,  a  CTD  was  collected  in  the  vicinity  of  the  R/P  FLIP  on 
Julian  Day  267,  up  to  a  depth  of  3885  m.  Another  CTD,  collected  at  37°04.19'  N 
134° 46.75'  W  on  Julian  Day  264,  was  used  to  extend  the  FLIP  CTD  down  to  the  nominal 
bottom  depth  at  the  R/P  FLIP,  4667  m.  The  sound  speed  profile,  based  on  the  CTD,  was 
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derived  using  the  UNESCO  equations  relating  conductivity  to  practical  salinity  and 
sound  speed  to  temperature,  salinity  and  pressure  [Fofonoffl983j. 

XBT  casts  were  collected  at  the  R/P  FLIP  before  and  after  the  fixed  station 
transmissions  on  Julian  Day  270  at  13:00  GMT  (a  400  m  XBT)  and  on  Julian  Day  272  at 
03:00  GMT  (a  750  m  XBT).  These  XBT  data  were  used  in  the  state  equation  of  [Mack¬ 
enzie  1981]  to  produce  sound  speed  profiles. 

This  information  allows  the  synthesis  of  the  sound  speed  profile  used  with  the 
ATLAS  model  to  produce  the  replica  vectors.  The  model  sound  speed  profile  is  plotted 
in  Figure  1  for  the  whole  water  column  and  for  the  top  750  m.  The  two  XBT  sound 
speed  profiles  are  also  plotted  in  dashed  lines.  A  30  m  deep  mixed  layer  was  observed  at 
the  R/P  FLIP,  followed  by  a  shallow  duct  above  100  m,  which  evolves  with  time,  as 
shown  by  the  two  XBT  sound  speed  profiles  taken  only  1 4  hours  apart.  The  deep  sound 
axis  as  at  a  depth  of  about  550  m. 

3.2,2.  Modeling  with  ATLAS 

The  replica  vectors  produced  by  the  AlLAS  normal  mode  model  are  based  on 
the  sound  speed  profile  in  Figure  1  and  a  bottom  loss  table  which  corresponds  to  a  GSM 
Bottom  Province  Type  3,  given  in  Table  I.  A  range  independent  medium  with  a  water 
depth  of  4667  m  is  considered  here. 

The  source  is  at  a  nominal  range  of  165.5  km  and  modeling  indicates  that  the 
VLA  is  in  the  trailing  edge  of  the  third  convergence  zone,  which  becomes  sharper  as  the 
source  is  raised  in  the  water  column  [Tranl989].  The  angular  spectrum  of  the  pressure 
field  produced  by  the  ATLAS  model  was  found  in  good  agreement  with  the  beamformed 
experimental  data.  The  replica  vectors  are  calculated  for  ranges  going  from  152.5  km  to 
177.5  km,  every  250  m,  for  sources  every  5  m  from  5  to  600  m. 
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Angle  (deg) 

Loss  (dB) 

0° 

0.0 

10° 

0.4 

20° 

3.2 

30° 

5.2 

40° 

6.8 

O 

O 

7.6 

60° 

8.2 

70° 

8.2 

80° 

8.2 

90° 

8.4 

Table  1:  Bottom  Loss  Table  in  dB  (GSM  Province  Type  3) 


Figure  I:  Fixed  station  sound  speed  profile. 
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3.3.  Matched  Field  Processing  Results 

Matched  field  processing  is  performed  with  the  Bartlett  processor  and  the  MV 
processor  with  stabilization  factors  equal  to  l(T\  i(T2  and  10'1.  Table  2  references  the  4 
minute  and  22  second  long  data  segments  used  to  estimate  the  covariance  matrix  for  each 
station.  The  sequence  of  three  source  depths  (300.  150  and  20  m)  provides  an  opportunity 
to  test  on  real  data  the  capabilities  of  matched  field  processing  to  localize  in  depth. 

For  the  three  source  depths,  the  covariance  matrix  is  computed  by  keeping  the 
frequency  bin  #  1844  in  the  2048-point  FFT  (200.195  Hz)  where  the  signal  is  the  strong¬ 
est  across  the  array. 


Tape  number 

GMT  data  segment 
start  time 

Source 
depth  (m) 

915 

20 

54:37 

300 

917 

21 

34:59 

150 

9TE 

21 

02:49 

20 

Table  2:  GMT  time  reference  for  the  covariance  matrix  estimate 


The  results  are  gray-level  displayed  in  Figure  2  to  7.  In  each  case,  the  gray  levels 
span  the  full  dynamic  range  of  the  ambiguity  surface  and  correspond  to  a  scale  of  dB  re 
upa.  Each  ambiguity  surface  is  normalized  by  the  level  of  its  maximum  peak  which  is 
characterized  by  its  calibrated  power  level  in  dB  re  uPa.  its  range  and  depth.  This  highest 
peak  can  be  used  to  localize  the  source.  Table  3  summarizes  the  results  of  matched  field 
processing  on  the  experimental  data.  First,  one  observes  that  the  ambiguity  surfaces  exhi¬ 
bit  their  maximum  peaks  near  the  nominal  source  range.  165.5  km.  One  notes  the  great 
granularity  of  the  range-depth  ambiguity  surfaces  and  the  hirge  number  of  sidelobes.  The 
Bartlett  processor  has  the  most  sidelobes  but  has  also  the  largest  dynamic  range.  The  MV 
processor  suffers  large  loss  due  to  mismatch,  the  replica  are  imperfect.  An  increase  in  the 
fraction  of  noise  y  for  the  MV  processor  results  in  more  robustness  [Cox  1988],  but 
reduces  the  relative  signal  level  in  the  ambiguity  surfaces.  A  stabilization  factor  of  10':  is 
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best  based  on  a  subjective  inspection  of  the  range-depth  ambiguity  surfaces. 

Wnile  localization  in  range  is  successful,  localization  in  depth  produces  mixed 
results.  When  the  source  is  at  300  m  depth,  the  Bartlett  processor  produces  a  highest 
level  peak  at  a  depth  of  195  m,  while  the  MV  processors  find  shallower  source  depths. 
When  the  source  is  at  a  depth  of  100  m,  the  Bartlett  processor  and  the  MV  processor  with 
a  stabilization  factor  y  of  l(T3  indicate  a  very  shallow  source.  On  the  other  hand,  the  sta¬ 
bilized  MV  processor  with  y  equal  to  1CT2  and  10'1  find  the  correct  source  depth  at  1 10  m. 
Finally  when  the  source  is  at  20  m,  all  the  processors  correctly  estimate  the  source  depth. 
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Table  3:  Summary  of  matched  field  processing  (fix^d  station  data) 
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Figure  4:  Matched  field  processing  on  the  data  of  Tape  917,  part  1 
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Figure  6:  Matched  field  processing  on  the  data  of  Tape  918,  part  1 . 
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3.4.  Matched  Field  Processing  Simulations 

In  an  effort  to  understand  the  results  of  the  matched  field  processing  on  the  real 
data,  simulation  results  are  now  presented  in  the  case  where  there  is  no  mismatch.  Then, 
causes  for  mismatch  are  analyzed  and  simulations  involving  array  tilt  are  performed. 

3.4.1.  No  Mismatch  Simulations 

A  perfect  situation  is  now  assumed  where  there  is  no  noise  and  no  mismatch 
(that  is  the  replica  vectors  computed  by  the  ATLAS  model  constitute  the  "truth").  A 
200  Hz  source  is  simulated  at  a  range  of  165  km  and  a  depth  of  300,  100  and  20  m.  The 
ambiguity  surfaces  for  the  Bartlett  processor  and  the  three  different  implementations  of 
the  MV  processor  are  computed  and  plotted  in  Figures  8  to  13. 

Since  there  is  no  mismatch  and  no  noise,  the  source  is  correctly  localized  in  all 
cases  and  the  dynamic  range  of  the  ambiguity  surfaces  is  much  larger  than  for  the  real 
data.  It  is  on  the  order  of  50  dB  for  the  Bartlett  processor  and  the  MV  processor  with  a 
stabilization  factor  y  of  l(T3,  and  on  the  order  of  40  and  30  dB  for  the  MV  processor  with 
stabilization  factors  equal  to  l(T2  and  10"1.  In  the  case  of  the  real  data,  there  is  loss  due  to 
mismatch  or  an  incorrect  modeling  of  the  wavefield.  The  Bartlett  processor  produces,  in 
all  these  simulation  cases,  a  large  number  of  sidelobes.  It  exhibits  the  usual  hat  shape 
also  observed  in  the  real  data  ambiguity  surfaces.  The  MV  processor  is  a  high  resolution 
processor  and  produces  a  single  peak  in  the  range-depth  cell  of  the  simulated  source. 
These  simulation  results  for  the  MV  processor  are  far  apart  from  what  is  observed  in  the 
real  data  case.  It  is  well  known  that  the  MV  processor  is  very  sensitive  to  mismatch  and 
lacks  robustness.  These  simulations  indicates  that  the  errors  in  the  replica  vectors  are 
sufficient  to  dramatically  perturb  the  MV  results. 

3.4.2.  Analysis  of  the  Mismatch 

There  are  several  potential  causes  for  the  large  mismatch  observed  in  the 
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matched  field  processing  runs  on  the  data.  Although  sound  speed  information  at  the 
R/P  FLIP  is  very  good,  one  does  assume  a  range  independent  medium.  This  may  not  be 
accurate  for  large  ranges  on  the  order  of  165  km.  There  is  no  available  range  dependent 
information  to  improve  the  processing.  Since  the  experiment  was  conducted  in  deep 
water,  the  bottom  should  have  only  a  weak  influence  in  the  processing  results.  Another 
cause  for  mismatch  is  the  position  errors  of  the  array  sensors.  The  VLA  array  was 
deployed  in  a  three  transponder  network  (long  baseline)  and  one  channel  per  array  sec¬ 
tion  was  navigated  during  the  experiment  [Sotirinl989].  Navigation  data  are  available  for 
Tape  915  and  917.  They  are  summarized  in  Figure  14  for  Tape  915  and  in  Figure  15  for 
Tape  917.  The  array  shapes  are  plotted  in  the  North-South  and  East-West  vertical  planes. 
The  navigation  data  indicate  that  the  array  behaves,  in  the  first  order,  like  a  pendulum 
and  that  the  largest  horizontal  displacement  across  the  array  aperture  is  bounded  by  15  m. 
Such  displacement  corresponds  to  a  1°  tilt.  Since  15  m  is  equal  to  two  wavelength  at  200 
Hz,  array  tilt  is  likely  to  have  a  major  impact  on  the  matched  field  processing  results.  It 
is  investigated  in  the  following. 
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Figure  10:  Matched  field  processing  simulation  (no  noise,  no  mismatch)  for  a  source 
depth  of  100  m,  part  1. 


Figure  11:  Matched  field  processing  simulation  (no  noise,  no  mismatch)  for  a  source 
depth  of  100  m,  pan  2. 


Figure  13:  Matched  field  processing  simulation  (no  noise,  no  mismatch)  for  a  source 
depth  of  20  m,  port  2. 
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3.4.3.  Tilted  Array  Simulation 

A  tilted  vertical  line  array  is  simulated  with  a  horizontal  displacement  of  15  m. 
The  source  is  assumed  at  165  km  from  the  array  top  element  and  at  a  depth  of  300,  100 
and  then  20  m.  The  tilted  array  simulation  geometry  is  summarized  in  Figure  16.  The 
replica  vectors,  which  were  used  earlier,  assume  a  straight  line  vertical  array,  thus  the 
impact  of  mismatch  due  to  array  tilt  is  assessed. 


a 

depth 


BOTTOM 


Figure  16:  Tilted  array  simulation 

The  ambiguity  surfaces  for  a  tilted  array  are  plotted  in  Figure  17  and  18  for  a 
source  at  300  m,  in  Figure  19  and  20  for  a  source  at  100  m,  and  in  Figure  21  and  22  for  a 
source  at  20  m.  The  matched  field  processing  results  are  summarized  in  Table  4.  Array 
tilt  results  in  mismatch  which  significantly  reduces  the  maximum  peak  power  and  shifts 
the  position  of  the  peak  in  range  and  depth  (the  probable  position  of  the  source).  The 
mismatch  also  reduces  the  dynamic  range  of  the  ambiguity  surfaces  to  a  few  dB  for  the 
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MV  processors.  The  estimated  position  of  the  source  is  off  in  range  by  4  to  7  km  from 
the  true  position  of  the  source.  The  300  m  deep  source  is  localized  at  a  depth  of  450  m 
while  the  100  m  deep  source  is  localized  at  a  depth  of  15  m.  These  simulated  ambiguity 
surfaces,  although  they  share  some  similarities  with  the  real  data  ambiguity  surfaces, 
correspond  to  a  worst  case  situation  where  array  tilt  is  maximum.  These  simulated  results 
indicate  that  array  tilt  can  be  held  reponsible  for  the  large  mismatch  observed  in  the  real 
data  ambiguity  surfaces. 
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Table  4:  Tilted  array  matched  field  processing  simulation 
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Figure  17:  Matched  field  processing  simulation  for  a  tilted  array  and  a  source  300  m 
depth,  part  1. 
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Figure  18:  Matched  field  processing  simulation  for  a  tilted  array  and  a  source  300  m 
depth,  part  2. 
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Figure  19:  Matched  field  processing  simulation  for  a  tilted  array  and  a  source  100  m 
depth,  parti. 
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4.  Source  Tow  Matched  Field  Processing 

4.1,  Introduction 

On  Julian  Day  272,  the  USNS  NARRAGAN3ETT  began  a  source  tow  for  10 
hours,  starting  at  00:00  GMT.  The  HX90  source  was  towed  at  a  nominal  speed  of  5  kts 
with  300  m  of  wire  out  (the  probable  source  depth  is  100  m),  and  was  projecting  a 
200  Hz  CVV  tone.  The  ship  was  tracked  by  radar  from  the  R/P  FLIP  during  the  first  five 
hours  of  the  ttansmission.  Data  were  recorded  continuously  on  25  tapes  from  Tape  966 
to  990  (each  tape  corresponds  to  23  minutes  and  28  seconds).  The  swell  heights  recorded 
at  the  R/P  FLIP  were  just  below  5  feet  and  wind  speeds  were  between  10  and  15  kts.  An 
XBT  cast  was  taken  on  the  R/P  FLIP  during  the  source  tow. 

4.2.  Source  Tow  Modeling 
4.2.1.  Sound  Speed  Profde 

The  data  from  the  XBT  cast  collected  on  the  R/P  FLIP  at  03:00  GMT  were  util¬ 
ized  with  the  state  equation  of  [Mackenziel981]  to  produce  the  sound  speed  profile  plot¬ 
ted  in  Figure  23.  The  environment  at  the  R/P  riMP  is  one  or  a  uoumc  duct  with  a  large 
and  well  defined  surface  duct  in  the  upper  hundred  meters.  The  estimated  depth  of  the 
towed  source  is  100  m,  right  at  the  boundary  of  the  two  waveguides.  Therefore,  this 
situation  may  be  difficult  to  model.  The  sound  speed  profile  across  the  water  column  is 
synthesized  by  using  a  desampled  version  of  the  750  m  XBT  sound  speed  profile  and  by 
blending  it  in  a  desampled  version  of  the  deep  CTD  at  the  R/P  FLIP  (the  one  used  in  the 
fixed  station  modeling).  The  desampled  XBT  sound  speed  profile  in  the  top  750  m  and 
the  synthetic  sound  speed  across  the  whole  water  column  are  also  plotted  in  Figure  23. 


Figure  23:  Source  tow  sound  speed  profile. 
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4.2.2.  Modeling  with  ATLAS 

As  in  the  case  of  the  fixed  station,  the  ATLAS  normal  mode  model  is  used  to 
produce  the  replica  vectors  for  the  source  tow  matched  field  processing  runs.  The 
environmental  parameters  are  similar  to  the  fixed  station  modeling  ones,  excepth  the 
sound  speed  profile.  A  range  independent  medium  is  assumed  and  a  bottom  loss  table 
corresponding  to  a  GSM  Bottom  Province  Type  3  is  used  (Table  1).  The  ATLAS  model 
predicts  the  first  convergence  zone  at  about  50  km.  A  fair  agreement  was  found  between 
the  arrival  structure  predicted  by  the  ATLAS  model  for  a  source  depth  of  100  m  and  the 
beamformed  data  [Tran  1989]. 

4.3.  Selection  of  the  Matched  Field  Processing  Stations 

4.3.1.  USNS  NARRAGANSETT  Tracking 

The  source  tow  began  at  00:00  GMT  on  Julian  Day  272.  During  the  first  four 
hours  and  fourty  five  minutes,  the  USNS  NARRAGANSETT  was  tracked  by  the  R/P 
FLIP  radar.  The  tracking  information  is  summarized  in  Figure  24.  The  range  of  the  ship 
as  a  function  of  time  is  derived  from  the  radar  tracking  information  and  plotted  in  Figure 
25.  The  dotted  line  in  Figure  25  allows  the  interpolation  of  the  ship  range  after  the  end  of 
the  radar  tracking. 
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Figure  24:  Radar  tracking  information. 
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Figure  25:  Ship  range  as  a  function  of  time 

4.3.2.  Selection  of  the  Stations 

Four  stations  are  processed  using  the  matched  field  processing  techniques.  They 
are  chosen  to  be  the  beginning  of  the  first  convergence  zone,  the  middle  of  the  first  con¬ 
vergence  zone,  the  trailing  edge  of  the  first  convergence  zone,  and  the  beginning  of  sha¬ 
dow  zone  folio  ving  the  first  convergence  zone.  The  stations  are  such  that  there  is  a  10 
km  range  increment  from  Station  #1  to  Station  #3,  and  from  Station  #3  to  Station  #4,  and 
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a  5  km  increment  from  Station  #1  to  Station  #2,  and  from  Station  #2  to  Station  #3.  Table 
5  summarizes  the  time  at  which  the  4  minutes  and  22  seconds  long  data  segment  begins 
for  each  station.  The  GMT  stall  time  of  the  data  segment  allows  the  estimation  of  the 
source  range  with  respect  to  the  VLA. 


Station 

Tape 

# 

GMT  Data 
Segment  Time 

Expected 
Range  (km) 

1 

980 

05:33:32 

40  to  45 

2 

981 

06:10:06 

45  to  50 

3 

983 

06:48:17 

50  to  55 

4 

986 

07:51:12 

60  to  65 

Table  5:  Source  tow  station  information 


4.4.  Doppler  Shift 

The  source  is  moving  away  from  the  VLA  receivers  at  a  nominal  speed  of  5  kts. 
The  VLA  array  is  essentially  fixed  since  deployed  from  the  R/P  FLIP  in  a  three  point 
moor.  The  Doppler  shifted  received  frequency  fr  is  given  by  [Campl970,  p  216] 


where  /,  the  transmitted  frequency  (i.e.  200  Hz),  c  the  sound  speed  and  v  the  tow  ship 
speed.  Assuming  a  constant  sound  speed  of  1500  m/s  and  a  tow  speed  of  5  kts,  the 
received  frequency  is  199.657  Hz.  It  is  different  from  200  Hz  by  1.4  times  the  bin  width 
of  the  Fourier  analysis  performed,  where  the  bin  width  is  0.244  Hz  (a  2048-point  FFT  is 
used  for  a  -250  Hz  to  250  Hz  frequency  band).  The  200  Hz  signals  from  the  towed  source 
is  strongest  in  bin  1843. 

Another  issue  is  the  differential  Doppler  shift  that  the  different  arrivals  can 
experience.  If  the  arrivals  have  a  relative  Doppler  shift  greater  than  the  FFT  binwidth, 
matched  field  processing  will  fail  to  coherently  recombine  the  multipaths.  The  differen¬ 
tial  Doppler  shift  is  easily  obtained  from  Equation  4.1  by  assuming  that  the  source  pro¬ 
ject  sound  with  a  given  take-off  angle  e  with  respect  to  the  horizontal,  while  the  source 
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has  only  a  horizontal  displacement.  While  there  is  no  Doppler  shift  for  sound  propagat¬ 
ing  vertically  (9  =  90°),  the  Doppler  shift  is  maximal  for  horizontally  propagating  sound 
(0  =  0°).  If  v  denotes  the  ship  horizontal  speed,  0  the  take-off  angle  of  the  sound  at  the 
source,  the  Doppler  shift  is  given  by 

A /  =  -  CQ--6--  (4.2) 

v 

cos©  +  — 
c 

v  is  replaced  in  Equation  4.1  by  its  projection  on  the  direction  0  (that  is  by  v  cos©)  to  pro¬ 
duce  Equation  4.2.  The  Doppler  shift  as  a  function  of  take-off  angle  is  plotted  in  Figure 
26.  The  Differential  shift  is  within  half  the  0.244  Hz  binwidth  for  all  practical  take-off 
angles  0  (i.e.  such  that  0  <70°). 


2 


Figure  26:  Doppler  shift  as  a  function  of  take-off  angle 

The  Doppler  shift  is  neglected  in  this  study  and  the  effects  due  to  frequency  vari¬ 
ations  are  not  considered  in  the  calculation  of  the  replica  vectors.  The  range  variations  of 
the  source  are  not  considered  either.  This  simple  approach  is  justified  in  the  following. 
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The  pressure  field  computed  by  the  ATLAS  normal  mode  model  can  be  written  as  [Gor¬ 
don  1984,  p  4] 


P  =  2  TT 

m  =  1 


2  n 
km  r 


Um(zs)  Um(z)e' 


(4.3) 


where  Nm  is  a  generalized  normalization  factor,  Um  the  m,h  mode  eigenfunction,  z  the 
depth  of  the  receiver,  zs  the  depth  of  the  source,  Km  the  horizontal  wavenumber  of  the  m'h 
mode  and  r  the  horizontal  range  from  the  source  to  the  receiver.  The  contribution  of  the 
mlh  mode  has  a  phase  given  by  Km  r  The  normalized  modal  eigenfunctions  are  gray-level 
displayed  in  Figure  27  for  the  whole  water  column  and  in  Figure  28  for  the  top  500  m.  A 
100  m  deep  source  excite  a  group  of  modes  around  mode  number  71  and  the  m‘h  mode 
discussed  in  the  following  can  be  assumed  to  belong  to  this  group  of  excited  modes.  For 
a  moving  source,  the  phase  of  the  m'h  mode  can  be  expressed  as 

<D  =  (Km  +dKm)(r  +  dr)  (4.4) 

where  dr  is  the  distance  over  which  the  source  moves  and  dKm  the  error  in  wavenumber 

due  to  a  slight  Doppler  frequency  shift.  The  phase  errors  when  using  kmr  can  be 

evaluated  by  keeping  only  the  first  order  terms 


<D-Kmr  =  Km  dr  +dKm  r 

.dr  . 

AO  =  Km  r( —  + - ) 

r  Km 

AO  _  dr  dKm 
O  ~  r  Km 


(4.5a) 

(4.5b) 

(4.5c) 


Since  the  source  moves  over  750  m  over  the  time  necessary  to  estimate  the  array  covari¬ 
ance  matrix  between  40  and  60  km,  the  relative  error  due  to  source  motion  in  range  is  at 
most  on  the  order  of  1.9  %. 

The  relative  phase  error  due  to  wavenumber  variations  is  shown  negligible  in  the 
following.  The  error  in  wavenumber  is 


dKm  =  Km(ti)  +  rfw)  -  Km(C0) 

dKm  =  — —  d  or 
d  co 


(4.6a) 

(4.6b) 
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where  co  is  the  circular  frequency.  Since  is  the  inverse  of  the  group  velocity  vg. 
Equation  4.6b  can  be  expressed  as 


The  Doppler  shift  d  co  is  given  by 


dKm  = 


d  co 


(4.7) 


where  V  is  the  speed  of  the  source,  cp  the  phase  velocity  of  the  m'h  mode, 
much  larger  than  V,  Equation  4.8  can  be  expressed  as  [Ozardl989] 


(4.8) 
Since  cp  is 


so  that 


d  k„ 


d  =- 


coV 


(4.9a) 

(4.9b) 


d  K_ 


(4.10) 


Here,  V  #  3  m/s  and  is  approximately  constant.  The  relative  error  in  horizontal 
wavenumber  computed  using  the  ATLAS  model  group  velocities  is  plotted  in  Figure  29 
against  mode  number.  It  is  bounded  by  0.18  %  and  an  order  of  magnitude  lower  than  the 
phase  errors  due  to  range  effects. 


This  preliminary  remarks  indicates  that  Doppler  shift  results  in  negligible  phase 
eirors  compared  to  source  range  variations.  Neglecting  Doppler  effects  in  the  replica  vec¬ 
tors  calculations  appear  justified. 


4.5.  Source  Tow  Matched  Field  Processing  Result 

Matched  field  processing  is  performed  using  the  Bartlett  processor  and  the  MV 
processor  with  stabilization  factors  equal  to  10‘\  10~2  and  KT1.  The  sequence  of  the  four 
selected  stations  provides  an  opportunity  to  test  on  real  data  the  capabilities  of  matched 
field  processing  to  localize  a  moving  source.  The  ambiguity  surface  are  computed 
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between  35  and  60  km  for  the  first  three  stations  and  between  47.5  and  72.5  km  for  the 
last  station.  As  before,  a  250  m  range  and  5  m  depth  cell  is  used.  The  maximum  depth  of 
the  ambiguity  surface  is  600  m.  The  results  are  gray-level  displayed  in  Figure  30  to  37. 
The  gray  levels  span  the  full  dynamic  range  of  each  ambiguity  surface  and  correspond  to 
a  scale  in  dB  re  pPa.  The  surfaces  are  normalized  by  their  highest  peak  power  which  can 
be  used  for  localization.  Table  6  summarizes  the  results  of  the  matched  field  processing 
on  the  source  tow  data. 

The  localization  in  range  is  excellent  since  the  ranges  of  the  highest  peak  in  the 
ambiguity  surfaces  correspond  to  the  predicted  ranges  for  each  station  given  in  Table  5. 
The  localization  in  depth  gives,  as  in  the  case  of  the  fixed  stations,  mixed  results.  The 
depth  of  the  highest  peak  is  either  shallower  or  deeper  than  the  expected  source  depth 
around  100  m.  Because  of  the  multiplicity  of  the  sidelobes,  it  is  difficult  to  estimate  pre¬ 
cisely  where  the  source  is.  The  dynamic  range  in  the  ambiguity  surfaces  indicates  a  large 
loss  due  to  mismatch  for  the  MV  processor. 
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Tape 

Bartlett 

Processor 

MV 

y  =  10‘3 

MV 

y  =  10"2 

MV 

y=  10’1 

980 

Depth  of  Max. 

200  m 

5Cm 

5'J  rn 

50  m 

Range  of  Max. 

41.5  km 

42.75  km 

42.75  km 

42.75  km 

Surface  Max. 

88.2  dB 

54.9  dB 

61.2  dB 

69.1  dB 

Dynamic  Range 

35  dB 

8.5  dB 

6.5  dB 

5  dB 

Figure 

30 

30 

31 

31 

981  . 

Depth  of  Max. 

50  m 

40m 

40  m 

40m 

Range  of  Max. 

44.25  km 

46.25  km 

46  km 

46  km 

Surface  Max. 

80.4  dB 

49.5  dB 

54.9  dB 

61.5  dB 

Dynamic  Range 

30  dB 

9  dB 

7  dB 

5  dB 

Figure 

32 

32 

33 

33 

983 

Depth  of  Max. 

135  m 

135  m 

20  m 

135  m 

Range  of  Max. 

51.75  km 

51.75  km 

50.75  km 

50.75  km 

Surface  Max. 

82.2  dB 

49.2  dB 

56.7  dB 

63.2  dB 

Dynamic  Range 

30  dB 

9  dB 

7  dB 

6  dB 

Figure 

34 

34 

35 

35 

986 

Depth  of  Max. 

15  m 

10m 

10  m 

10  m 

Range  of  Max. 

60.25  km 

60  km 

60.25  km 

60.25  km 

Surface  Max. 

78.2  dB 

49.2  dB 

54.8  dB 

62.0  dB 

Dynamic  Rang 

25  dB 

9.5  dB 

8  dB 

7  dB 

Figure 

36 

36 

37 

37 

Table  6:  Summary  of  the  matched  field  processing  (source  tow  data) 


4.6.  Matched  Field  Processing  Simulation 

As  in  the  case  of  the  fixed  station  processing,  simulation  results  are  now 
presented  to  help  the  interpretation  of  the  results  of  the  matched  field  processing  on  the 
real  data. 

4.6.1.  Analysis  of  Mismatch 

As  pointed  out  in  the  previous  section,  the  results  show  a  somewhat  large 
mismatch.  The  source  tow  matched  field  processing  runs  are  performed  on  a  much 
shorter  range  than  the  fixed  stations,  therefore  the  assumption  of  a  range  independent 
medium  appears  justified.  The  sound  speed  profile  is  derived  from  an  XBT  cast  taken 
during  the  source  tow  and  sound  speed  mismatch  is  unlikely.  The  impact  of  the  source 
motion  on  the  ambiguity  surfaces  is  investigated.  Like  in  the  case  of  the  fixed  stations, 
array  tilt  is  expected  to  have  a  major  impact  on  the  performance  of  the  matched  field 
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processing  techniques  and  will  also  be  discussed. 


4. 6.2.  Simulation  of  the  Tow 


A  source  tow  simulation  for  Station  #2  is  performed,  where  the  source  is 
assumed  at  a  depth  of  100  m  in  a  noise-free  environment  and  moves  away  from  the  VLA 
array.  4  minutes  and  22  seconds  worth  of  data  are  used  to  estimate  the  real  data  covari¬ 
ance  matrix  thus  for  a  ship  speed  of  5  kts  or  2.75  m/s,  the  source  displacement  is  721  m 
or  3  range  cells  (each  cell  is  250  m  wide).  The  100  m  deep  moving  source  is  simulated 
by  extracting  the  unnormalized  complex  ATLAS  pressure  fields  that  correspond  to 
ranges  equal  to  47.5  km,  47.75  km  and  48  km.  Denoting  those  complex  vectors  by  X0,  Xl5 
X2,  the  covariance  matrix  is  computed  by  using 


R  = 


3  Xq  X0 


X0 X0"  +  X,  X"+X2  X2W 


The  replica  vector  that  corresponds  to  X,  is  given  by 


X, 


(X,WX,)‘- 


(4.11) 


In  this  simulation,  the 


"true"  source  position  is  considered  at  47.5  km  from  the  VLA.  Matched  field  processing 
is  performed  on  this  covariance  matrix  and  the  ambiguity  surfaces  are  plotted  in  Figure 
38  and  39.  The  only  mismatch  is  due  to  the  averaged  covariance  matrix.  The  ambiguity 
surfaces  are  similar  to  the  ones  of  Figure  40,  which  correspond  to  the  perfect  case  where 
there  is  no  mismatch  (a  non  moving  source  is  assumed  at  47.5  km  and  100  m  depth). 
The  highest  peaks  in  Figure  38  and  39  are  at  the  right  range-depth  cell  but  their  power  is 
reduced  by  almost  3  dB  for  the  Bartlett  processor  and  almost  5  dB  for  the  MV  processor. 
The  pressure  field  sampled  by  the  VLA  is  not  stationary  since  the  source  moves  away  in 
range.  This  simulation  shows  that  averaging  over  time  snapshots  to  estimate  the  covari¬ 
ance  matrix  results  in  a  reduction  of  the  maximum  peak  power  and  a  slight  broadening  of 
the  main  peaks.  But  it  also  shows  that  averaging  does  not  create  enough  mismatch  to  des¬ 
troy  the  ambiguity  surfaces  and  to  reduce,  in  a  dramatic  way,  their  dynamic  range.  Thus, 
the  amount  of  averaging  proposed  in  Section  2  and  used  to  estimate  the  covariance 
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matrix  appears  acceptable. 

To  illustrate  the  undesirable  effects  of  using  a  larger  number  of  averages  (i.e. 
using  a  longer  data  segment)  to  estimate  the  array  covariance  matrix,  the  ambiguity  sur¬ 
faces  for  a  source  moving  over  1500  m  or  6  range  cells  are  calculated  and  plotted  in  Fig¬ 
ure  41  and  42  This  simulates  the  case  where  252  dyadic  vector  outer-products  are  used  to 
compute  the  covariance  matrix  with  a  2048-point  FFT  operating  on  the  data.  These 
ambiguity  surfaces  show  a  large  peak  power  reduction  due  to  mismatch,  almost  5  dB  for 
the  Bartlett  processor  and  almost  8  dB  for  the  MV  processors.  Furthermore  the  ambiguity 
surface  for  the  Bartlett  processor  appears  qualitatively  quite  different  from  the  one  with 
no  mismatch. 

4.6.3.  Simulation  of  Tilt 

As  in  the  case  of  the  fixed  stations,  the  array  was  navigated  in  the  three  tran¬ 
sponder  net.  The  shapes  of  the  array  in  the  North-East  and  East-West  vertical  planes  as 
well  as  from  above  in  a  horizontal  plane  are  plotted  in  Figure  43,  44  and  45  for  Tapes 
980,  983  and  986.  The  navigation  shows  that  the  array  is  tilted  with  maximum  horizontal 
displacement  from  top  to  bottom  on  the  order  of  15  m.  The  array  was  fairly  straight  when 
Tape  986  was  recorded  and  mismatch  due  to  tilt  should  minimal  in  this  case. 

A  tilted  array  is  simulated  at  a  range  of  47.5  km  from  the  VLA  in  Figure  46  and 
47.  The  source  is  assumed  at  a  depth  of  100  m.  There  is  no  noise  and  no  other  mismatch 
than  array  tilt.  As  expected,  there  is  a  large  loss  due  to  mismatch  for  the  MV  processors. 
The  power  of  main  peak  is  significantly  reduced  as  well  as  the  dynamic  range  in  the 
ambiguity  surface.  The  localization  in  range  is  underestimated  by  1.75  km  and  the 
source  depth  is  estimated  at  95  m  instead  of  100  m. 

4.6.4.  Tow  and  Tilt  Simulation 

Since  averaging  over  space  and  array  tilting  occur  simultaneously  in  the  case  of 
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the  real  data,  a  simulation  involving  both  phenomenon  is  presented  for  all  four  stations. 
A  tilted  array  is  simulated  with  a  15  m  horizontal  displacement  from  top  to  bottom,  with 
the  source  at  100  m  depth  and  moving  over  3  range  cells  starting  at  42  km  for  Station  #1 , 
at  47.5  km  for  Station  #2,  at  52  km  for  Station  #3  and  at  62  km  for  Station  #4.  The 
covariance  matrix  is  computed  by 


R  = 


2  xr.x« 

i  =  1 

x/fxr; 


(4.12) 


where  r,  belongs  to  the  set  of  ranges  {42,  47.5,  52,  62}  and  r,  =rj+0.25  (/ — 1),  i*  1.  The 
ambiguity  surfaces  are  computed  in  each  case  for  the  Bartlett  processor,  and  the  three 
implementations  of  the  MV  processor.  Good  similarity  is  observed  with  the  experimen¬ 
tal  runs  for  the  first  three  stations.  For  Tape  986  or  Station  #4,  the  simulation  results  are 
different  from  the  experimental  run.  This  is  not  surprising  since  the  array  tilt,  in  this  case, 
is  minimal.  The  simulation  results  are  summarized  in  Table  8.  The  simulation  shows  that 
array  tilt  coupled  with  averaging  over  range  results  in  poor  estimate  of  th^  source  depth 
for  the  three  first  stations.  The  range  estimate  is  generally  close  to  the  true  range.  It  is 
generally  off  by  1.75  km  except  for  Station  #4.  The  matched  field  processing  performed 
on  Station  #4  produces  a  good  estimate  of  the  source  depth  but  a  poor  estimate  of  the 
source  range  (the  error  is  either  4.75  km  or  5  km).  Also,  the  MV  processors  suffer  the 
most  loss  due  to  mismatch  and  the  most  reduction  of  dynamic  range  in  the  ambiguity  sur¬ 


faces. 
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Station 

Processor 

Bartlett 

MV 

MV 

MV 

y=  10'3 

y-  10-2 

y  =  10“' 

1 

Depth  of  Max. 

70  m 

60  m 

60  m 

60  m 

Range  of  Max. 

41.25  km 

41.5  km 

41.5  km 

41.5  km 

Surface  Max. 

-5dB 

-47.8  dB 

-37.8  dB 

-27.8  dB 

Dynamic  Range 

40d3 

3  dB 

3  dB 

3  dB 

Figure 

48 

48 

49 

49 

2 

Depth  of  Max. 

95  m 

35  m 

35  m 

35  m 

Range  of  Max. 

46  km 

46.75  km 

46.75  km 

46.75  km 

Surface  Max. 

-4.95  dB 

-47  dB 

-37  dB 

-27  dB 

Dynamic  Range 

50  dB 

4  dB 

4  dB 

4  dB 

Figure 

50 

50 

51 

51 

3 

Depth  of  Max. 

55  m 

20  m 

20  m 

20  m 

Range  of  Max. 

50.75  km 

50.25  km 

50.25  km 

50.25  km 

Surface  Max. 

-6  dB 

-46.2  dB 

-36.2  dB 

-26.2  dB 

Dynamic  Range 

45  dB 

4.7  dB 

4.7  dB 

4.7  dB 

Figure 

52 

52 

53 

53 

4 

Depth  of  Max. 

1 10  m 

110  m 

110m 

110  m 

Range  of  Max. 

66.25  km 

67  km 

67  km 

67  km 

Surface  Max. 

-9  dB 

-50  dB 

-40  dB 

-30  dB 

Dynamic  Range 

38  dB 

1.1  dB 

1.1  dB 

1.1  dB 

Figure 

54 

54 

55 

55 

Table  8:  Summary  of  the  coupled  tilt  and  tow  simulation 


5.  Conclusions 

This  study  demonstrates,  using  real  data,  that  matched  field  processing  is  a  viable 
processing  technique  at  200  Hz  with  a  900  m  long  vertical  line  array.  Because  of  the  high 
source  level  used  in  this  transmission  experiment,  the  source  is  easily  detected  and  detec¬ 
tion  issues  are  not  discussed.  This  report  demonstrate  the  good  capabilities  of  range 
localization  of  a  source  in  deep  water  and  for  the  200  Hz  mid-frequency.  The  environ¬ 
mental  information  collected  at  the  time  of  the  experiment  allows  an  appropriate  model¬ 
ing  of  the  acoustic  transmissions  using  the  ATLAS  normal  mode  model.  While  the 
acoustic  modeling,  still,  is  imperfect,  the  source  can  be  localized  with  good  accuracy  in 
range.  Mismatch  exists  mainly  because  of  a  1°  tilt.  This  tilt  is  severe  at  the  operating  fre¬ 
quency  since  it  corresponds  to  twice  the  wavelength  at  200  Hz.  Simulations  show  that 
array  tilt  causes  slight  range  localization  errors  but  potentially  large  errors  for  depth 
localization.  Source  depth  estimation  appears  more  sensitive  to  mismatch. 
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Figure  32:  Matched  field  processing  on  the  data  of  Tape  981 ,  part  1 
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Figure  35:  Matched  field  processing  on  the  data  of  Tape  983,  part  2 
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Figure  36:  Matched  field  processing  on  the  data  of  Tape  986.  part  1 


47.5  50.0  52.5  55 .0  57.5  60.0  62.5  6: 


f'o .  0  7 


MV  Processor  Ronge-Depth 
Rmbiguity  Surface 

Surface  Max.  :  54.8  dB 

Range  of  Max.:  60.25  Un 
Depth  of  Max .  :  10  m 

Fraction  of  Noise:  0.01 


Range  (Un) 


-5  -4  -3 

Relative  Power  (d&) 


&  40 
a 


47.5  50.0  52.5  55.0  57.5  60.0  62.5  65.0  67.5  70.0  7 


MV  Processor  Ranee-Depth 
Rmbiguity  Surface 
Surface  Max.  :  62  dB 


Range  c;  Max. 
Depth  cf  Max . 


69.25  Un 
19  m 


■rcction  of  Noise:  0.  1 


Rcnpe  (Un) 


■7  -6  -5  -4 


Relative  Power  (d&> 


Figure  37:  Matched  field  processing  on  the  data  of  Tape  986,  part  2. 
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Figure  38:  Source  tow  matched  field  processing  simulation,  part  1 
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Figure  41:  Source  tow  matched  field  processing  simulation  (long  averaging  time),  part 
1. 
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Figure  42:  Source  tow  matched  field  processing  simulation  (long  averaging  time),  part 
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Figure  45:  Array  navigation  for  Tape  986 
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Figure  49:  Matched  field  processing  simulation  (tilted  array  and  towed  source  at  4 
km),  part  2. 
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Figure  50:  Matched  field  processing  simulation  (tilted  array  and  towed  source  at  47.5 
km),  part  1. 
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Figure  54:  Matched  field  processing  simulation  (tilted  array  and  towed  M>urce  at  62 
km),  pan  1. 
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